Astronomy & Astrophysics manuscript no. messineo 


©ESO 2011 


November 1, 2011 





Near- and Mid-Infrared colors of evolved stars in the Galactic 
plane. The Ql and Q2 parameters. 

M. Messineo' and K. M. Menten' and E. Churchwell^ and H. Habing^ 



' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

e-mail: me s s ineoOmpi f r -bonn .mpg.de 
- Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands 

' Department of Astronomy, University of Wisconsin-Madison, 475 N. Charter street, Madison, WI 53706, USA 
' Received September 15, 1996; accepted March 16, 1997 

(N ■ 

ABSTRACT 

o 

Context. Mass-loss from evolved stars chemically enriches the interstellar medium (ISM). Stellar winds from massive stars and their 
explosions as supemovae shape the ISM and trigger star formation. Studying evolved stars is fundamental for understanding galaxy 
formation and evolution, at any redshift. 
■ Aims. We aim to establish a photometric identification and classification scheme for Galactic mass-losing evolved stars (e.g., WR, 

RSG, and AGB stars) with the goal of identifying new ones, and subsequently to use the sample as tracers of Galactic structure. 
Methods. We searched for counterpai'ts of known Galactic WR, LBV, RSG, and 0-rich AGB stars in the 2MASS, GLIMPSE, and 
, MSX catalogs, and we analyzed their properties with near- and mid-infrared color-color diagrams. 

O Results. We used the Ql parameter, which is a measure of the deviation from the interstellar reddening vector in the J - H versus 
' H-Ks diagram, and we defined a new parameter, Q2, that is a measure of the deviation from the interstellar reddening vector in the 

(~| I J-Ks versus A^s-[8.0] diagram. The latter plane enables to distinguish between interstellar and circumstellar reddening, and to identify 

stars with circumstellar envelopes. WR stars and late-type mass-losing stars (AGBs and RSGs) are distributed in two different regions 
I . of the Ql versus A^s-[8.0] diagram. A sequence of increasing [3.6] - [4.5] and [3.6] - [8.0] colors with increasing pulsation amplitudes 

O ' (SRs, Miras, and OH/IR stars) is found. Spectra of Miras and OH/IR stars have stronger water absorption at 3.0 /jm than SR stars or 

most of the RSGs. Masing Miras stars have water, but stronger SiO (~ 4 um) and CO2 absorption (~ 4.25 /jm), as suggested by their 
C/3 , [3.6] - [4.5] colors, bluer than those of non masing stars. A fraction of RSGs (22%) have the bluest [3.6] - [4.5] colors, but small Q2 

values. We propose a new set of photometric criteria to distinguish among IR bright Galactic stars. 

Conclusions. The GLIMPSE catalog is a powerful tool for photometric classification of Galactic mass-losing evolved stars. Our new 
criteria will yield many new RSGs and WRs. 

^ ' Key words. Stars: mass-loss, dust: extinction. Stars: Wolf-Rayet, Stars: supergiants, Stars: late-type, Galaxy: stellar content 

00 

\^ 1 . Introduction ity of hot massive stars are not associated with dusty envelopes; 

WRs produce only 5% of the dust lGehrz,,1989.) . While mas- 
^ Evolved stars ai-e intrinsically bright and, when their luminosi- sive stars are typically f ound in or near HII regions a nd well 
^ ties can be established, they may serve as distance indicators, ti-ace galactic spiral ai-ms ("Georg elin & GeorgelinI 19761) . AGBs, 
I Thus, they are ideal tracers of Galactic structure. Moreover, which are older and dynamically relaxed, may serve as di- 
• . ! since they lose mass at high rates, they chemically enrich the rect ti-acers of the Galactic gravitational potential (iHabing et al.1 
^ ' Galactic interstellar medium (ISM). 2006) . Kinematics of AGBs have been successfull y used to map 

mass components devoid of gas (e.g, galactic bars lHabing et al.l 
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_ Intermediate mass stars (from ~ 2.2 to ~ 8 M^) lose mass at MlMessineo et ^l2002h . 
jlj , high rates when they go trough the Asymptotic Giant Branch 

phase (AGB stars), and, thereafter, the planetary nebula (PN) A census and classification of Galactic evolved stars is of 

phase, before cooling as white dwarfs. Evolved massive stars primary importance to understand galaxy formation and evolu- 

(> 8 Mq) go through several short-lived evolutionary phases tion. AGBs and RSGs are still poorly modeled, and their relative 

(blue supergiants, BSGs, luminous blue variables, LBVs, Wolf- contribution to the integrated light of a galaxy is quite uncer- 

Rayet stars, WRs, yellow supergiants, YSGs, and red super- tain. Major difficulties arise from our poor understanding of the 

giants, RSGs). Eventually, they explode as supernovae leaving complex process of mass-loss. Moreover, observational identifi- 

a neutron star or a black-hole. Late-type stars (red giants, AGBs, cation of evolved stars is a difficult task. Optical studies cannot 

and RSGs) ar e the most prolific dust producers in the Universe penetrate the inner Galactic regions due to interstellar extinc- 

dGehr j[T989l). The archetype VY CMa has a mass-loss rate of tion; stellar classification in the near-infrared is more difficult, 

lO-'* Mo y-i (iHarwit et al.il200Tl) . While RSGs may have higher because intrin sic near-infrare d colors of evolved stars span a 

mass-loss, since M cc L, luminous red giants and AGB are more small range ( Koornnee il983h . Spectroscopic observations may 

numerous and, account for more than 70% of the dust in the in- solve the ambiguity, but they requires more telescope time than 

terstellar medium. M and L are the stell ar mass loss rate and lu- photometry, and are feasible only for a limited number of ob- 

minosity, respectively (e.g. lGehr3ll989h . In contrast, the major- jects. Photometric variability, mass-loss, interstellar extinction, 



1 



Please give a shorter version witii: \authorrunning and/or \titilerunning prior to \maketitle 



and poor knowledge on stellar distances hamper the detection 
and classification even of these brightest stars. 

Several near-infrared and mid-infrared surveys of the whole 
Galactic plane are now (or will soon beco me) available, e.g . 
the Two Micron All Sky survey (2MASS) ("Cutri et al ' 
the D eep Near-Infrared Survey (DENIS) (Epchteinet al 
ffhe 2005), the UKIRT Infrared Deep Sky Survey (UKIDSS) 
(Lucas etal. 2008), the VISTA survey ( Minnitietal 20im , 
the Midcourse Space Experiment (MSX) jPrice et al.l 120011) . 
the ISO infrared su rvey of the Galactic Plane (ISOGAL) 
dSchuUeret al.ll2003h . and the Galactic Legacy Infrared Mid- 
Plane Survey Extraordinaire (GLIMPSE) ( Benjamin et alj2003l: 
IChurchwell et al. 2009). The combined use of near- and mid- 
infrared measurements permits detection of evolved mass-losing 
stars, as d emonstrated with ISOGAL and MSX data (e.g. 
lAIard et al.l 12001: Messineo et al. 2005). However, the classifi- 
cation of evolved stars using GLIMPSE data is still largely unex- 
plored. Simi lar studies have been carried out on the Magellanic 
Clouds (e.g. iBonanos et al.ll2010t [Buchanan et al.ll2006l l2009t 
lYang& Jiangl l201 l|). However, these studies suffer from small 
number statistics, and are limited to an environment with low 
metallicity. 

In the Milky Way, we know about 500 RSGs, 350 WR 
stars ([van der Hucht 2001), and a dozen confirmed LBVs 
(IClark et al.l i2005). Several thousands of AGBs have been de- 
tected via t heir maser emission, or photom etric pulsatio n prop- 



erties (e.g. lAlard et all 120011: iDeguchi et al. 2004: 



l200lHHabing et alJl2006l:lMessineo et anr2002:,Sevensteill2002 



and references therein). These numbers are too small for quan- 
titatively constraining theories of Galactic formation and e volu- 
tion (e.g. iHabinget all 120061: IVauterin & Deionghel[T998h . and 
for significantly sampling short-lived evolutionary phases and 
their contribution to Galactic chemical enrichment. Galactic 
models predict about 9 m illion Miras, ~ 5000 M supergiants, 
and ~ 3000 WR stars (e.g. lGehrdll989l) . 

Stellar classification based on infrared two-color diagrams 
was successfully established in the late 1980's, based on data 
from the Infrared Astronomical Satellite Point Source cat- 
alogue (the IRAS PSC) (e.g. van der Veen & Habing 1988). 
The 60//m/25//m versus 25/im/12/im color-color diagram distin- 
guishes between a carbon and an oxygen-rich chemistry-, and 
shows a sequence of O-rich circumstellar shells with increas- 
ing mass-loss rates. Successful selection of mass-losing late- 
type stars has b een made wi th data from the MSX survey (e.g. 
Messine o et al.l l2b04. 2005; ISevensteiil2002l: ISiouwerman et alj 
2009). A dust sequence of O-rich envelopes is seen also with 
MSX colors (Sevensterl [19991 ISiouwerman et al. 2009). MSX 
mapped the whole Galactic Plane within 5° of latitude (plus the 
IRAS gaps) in six bands (the Z?l-band is centered at 4.3 yum, 
the B2 band at 4.35 yum, the A-band at 8.28 yum, the C-band at 
12.13 yum, the D-bandat 14.65 yum, and the E-hand at 2 1.34 /^m). 
The A-ba nd has a sensitivi ty of 0.1 Jy, and a spatial resolution 
of 18.3" (lPriceetalJl200lh . We indicate with [A],[C],[D], and 
[E], the magnitudes in the 4 MSX bands at 8.28 fim, 12.13//m, 
14.65yum, and 21.34yum, respectively. 

The IRAC camera on board the Spitzer Space Telescope of- 
fers a new view of the sky in four filters (bands are centered 
at 3.6, 4.5, 5.8, and 8.0 fj.m iFazio et all l2004 . with a spa- 
tial resolution 6 times better than MSX, and a better sensitiv- 
ity. We indicate with [3.6], [4.5], [5.8], and [8.0], the magni- 
tudes in the four GLIMPSE filters. Buchanan et al.l (12009) and 
[Bonanos et al. (2010) developed a set of color criteria for clas- 
sifying luminous 8 fj.m sources in the Large Magellanic Cloud 
(LMC). HII regions can be isolated by their much redder K^^-A 



colors, while measurements in the four IR AC bands dist i nguish 
between RSGs, O-rich and C-rich AGBs. lYang & Jiand (1201 ih 
revised the pulsation properties of R SGs in the L MC using a 
comprehensive list of 191 RSGs (Buchanan et al. 2009). Only 
24% of RSGs show regular pulsation ( Yang & Jiang,.201 1). The 
IR selected sub-sample of LMC RSGs appears to have redder 
^rs-[8.0] colors (likely due to higher mass-loss rates) than the 
average sample. RSGs are found to have bluer [3.6] - [4.5] col- 
ors than other luminous red stars. 

So far, little work has been done on Galactic evolved stars 
based on GLIMPSE data. WRs have an inf rared excess, due t o 
free-free emission f rom dense stellar winds dCohen et al.lll975h . 
iHadfieW et all (12007) established a successful selection of WRs 
based on a combination of 2 MASS and GLIMPSE co lors, which 
has been recently revised by iMauerhan et al.l ( 1201 lb . WRs can 
be identified by combining three color-color diagrams, the J - 
H versus J-K, diagram, the [3.6] - [8.0] versus [3.6] - [4.5] 
diagram, and the ^rs-[8.0] versus J-K^ diagram (Figs. [1] [3] and 
IS. Known WRs have [3.6] - [8.0] > 0.5 mag and [3.6] - [4.5] > 
0.1 mag, and in the A's-[8.0] versus J-Ks diagrams are located 
on a different sequence than normal OB stars: their distribution 
with redder colors does not follow the reddening vector. Hadfield 
et al. reported a detection rate of 65% for emission line massive 
stars, and 7% for WR stars. Mauerhan et al. reached a detection 
rate of 95% for emission line massive stars, and up to 40% for 
WR stars using additional constraints based on the versus 
J-Ks diagram, and X-ray emission. 

An overview of near and mid-infrared colors of bright mass- 
losing stars in the Milky Way is missing, and providing one is the 
purpose of the present work. In this work, we study near-infrared 
and mid-infrared color properties of known Galactic evolved 
stars-, and define salient color-criteria for selecting new bonafide 
WR stars, RSGs, and AGBs. In Sect. 2, we report on previous 
analyses of evolved stars, which are based on GLIMPSE data. 
In Sect. 3, we describe the samples of known evolved stars, on 
which we base the new color-criteria. In Sect. 4, we describe the 
assumed interstellar extinction ratios. Color properties are ana- 
lyzed in Sect. 5, and discussed in Sect. 6. 

2. Sample selection 

We collected samples of known evolved massive stars (RSGs, 
WRs, and LBVs) with spectroscopically determined spectral 
types. The lists of AGBs were compiled on the basis of maser 
properties and/or photometric variability information (SiOmas- 
ing stars, OH/IR stars, Miras, and semiregular stars iHabingl 
1996). 

Throughout the text we use the terms of late-type stars and 
early-type stars. These definitions are only based on effective 
temperatures and not on luminosity classes. AGBs and RSGs 
are late-type stars, while LBVs and WRs are early-type stars. 

For every object, we searched for counterparts in the 2MASS 
All-Sky Catalog of Point Sources ( Cutri et al. 2003), in the third 
release of DENIS data av ailable at CDS (catalo g B/denis) and 
in the GLIMPSE catalog (ISpitzer Sciencdl2009t) using a search 
radius of 2". The 11/293 (GLIMPSE) catalog available from 
CDS merges the three surveys GLIMPSE-I (v2.0), GLIMPSE- 
II (v2.0), and GLIMPSE-3D: Catalog and Archive records are 
also merged. 

We searched for counterp arts in the Version 2.3 of the MSX 
Point Source Catalog (PSC) (lEgan et al.l l2003l) using a search 
radius of 5". 

Only sources detected in at least one mid-infrared band were 
retained for further analysis. Only 2MASS sources with good 
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photometry were retained (i.e., with red flags 1, or 2, or 3, and 
quahty flags A, B, C or D). Specific details on each sample are 
given in the following sub-sections. 



2.1. Wolf Rayet Stars 

Stars with initial masses larger than 40 Mq enter the WR phase 
when their H surface fractional mass abundance falls below 0.3 
and their surface temperature is roug hly above 20,000 K (e.g. 
iBressan et anil984l: iFiger et al.lll997h . WRs lose mass at high 
rates (~10"^ Mq y ') with their strong ionized winds, in which 
free-free emission is produced. This emiss ion can be detecte d 
as an infrared excess long-ward of 2 /im (ICohen etalJll975h . 
Their spectral energy distribution can be described by a cen- 
tral star plus an infrared excess, well described by a power law 
with a spectral index, a, of -0.6, which corresponds to the pre- 
dicted index for a free-free emission generated in stellar winds 
(iFelli & Panagial[T98TllWright & Barlow|[T97l . 

WC stars show circumstellar du st, that their d ust is likely due 
to interactions in multiple systems (IWatersll2010h . 

We searched for GLIMPSE and MSX counterparts of 
a sampl e of 345 Ga l actic WR stars (iMauerhan et al.l 201 it 
[Messineoetin 1201 U 120091; IShara et al.l 1201 U Ivan der Huchtl 
[200 1). Coordinates by Mauerhan et al. and Messineo et al. 
are from 2MASS, while those from Shara et al. are from the 
The Naval Observatory Merged Astrometric Dataset (NOMAD) 
dZach arias et al.l 12004*). The positions of WR stars listed 



[8.01-band emissi on of two LBVs ( e.g. iPasquali et all 120021: 
iPeeters et al.ll2002HUmana et alJl2010[) . 

We considered a list of 42 confirmed or can didate LBVs 
in flie Milkv NVav (IClark et all 12 005: Gvaramadze et alJl20ld : 



dZach i 
by B'a 



by [van der Hucht (2001) are from astrometric catalogs (e.g. 
Hipparcos), but for a few stars less accurate positions are given. 
Initially, we searched for possible associations within 10"; a few 
matches larger than 2" were individually checked , and coordi- 
nates updated by using the spectroscopic catalog bvl Skifil(l2010h . 
The GLIMPSE and 2MASS astrometries agree within 1". 

We retained the 120 GLIMPSE matches within 2", and the 
70 MSX matches within 5". The low number of GLIMPSE 
matches is mainly due to the survey coverage (longitude / < 65° 
and latitude b < 1°). 195 out of 345 WRs are located in the 
range (j/j < 65° and \b\ < 1°). 22 WRs stars have both MSX 
and GLIMPSE 8 yum measurements; the mean value of their 
([8.0] - [A]) color is -1-0.22 mag with a standard deviation of 0.5 
mag. Discrepant measurements, more than 2.5 sigma, are found 
for WR 125 and WR 93. Both sources are colliding wind bina- 
ries of the same type (WC7 + 09) (Norci et aL.200Z) . For the 
du sty WR125 WR, varia ble mid-infrared flux has been reported 
bv lWilliams et al.1 (11994 . 



2.2. Luminous blue variables 

LBVs are rare massive s tars in transition toward t he Wolf-Rayet 
phase (e.g. [C lark et al.l 120051: ICo nti et al. 19951 iMartins et all 
l2007t iNota et al.lll995FlSmifli et al.l i2004). They are character- 
ized by cyclic photometric variations with amplitudes of 1-2 
mag and yearly timescales, at nearly constant bolometric lu- 
minosity. Giant eruptions with changes in visual brightness by 
more than 3 mag may also occur LBVs have high mass-loss 
rates; the archetypical LBV, 77 Car, has a mass-loss rate of ~ 10"** 
Mq yr"' (e.g. Abraham et al. 2005). 

LBVs may be surrounded by extended dusty nebulae, which 
show a variety of morphology and composition. Amorphous and 
cristalline silicate, forst erite, as well as iron grains, have been 
observed ( Lamers et a I.I I2OOII: IPeeters et all 120021: lUmana et alj 
120 lot [Waters. ,2010) . PAH emission are found to dominate the 



IMauerhan et alJ2010trMessineo et aL201 L.2009.) . We found 18 
counterparts in the GLIMPSE catalog and 25 in the MSX cat- 
alog. Five stars have 8 yum measurements in both the MSX and 
GLIMPSE surveys. However, MSX measurements of WR102ka 
were di scarded because of h i gh crowding iri the g alactic center 
region (" Barniske et"an l2008t iHomeier et alJ 120031) . For the re- 
maining four sources, the mean value of their [8.0] - [A] color is 
-1-0.15 mag, with a standard deviation of +0.19 mag. 



2.3. Red supergiant stars 

RSG stars are late-type stars of ini tial masses between 8 and 40 
Mq, burning He in their core (e.g. lLevesque etalj|2005h . RSGs 
have luminosities log(L/L0) from 4 to 5.3, often show irregu- 
lar photometric variability, and are surrounded by dusty circum- 
stellar envelopes. T heir mass-loss rates range from ~ 10~^ to 
~ 10 ^^Mo yr-i (e.g. lJosselin et al.ll2000l: I Verhoelst et al.ll2009l) . 

Over 1000 stars are li s ted as RSGs in the spectroscopic 
catalog compiled by ISkifll 12010). However, several measure- 
ments are from the early 1900s, and need to be confirmed. 
A significant number of RSG stars were recently discov- 
ered as members of young massive clusters (e.g. 2MASS and 
GLIMPSE Messineo et al. 2009). Five young massive stellar 
clusters (RSGCl, RSGC2, RSGC3, RSGC4, and RSGC5), ex- 
traordinarily rich in RSGs (14, 26, > 8, > 9, 7), have been lo- 
cated between 25° and 30° of longitude, at a distanc e of about 6 



kpc probably at the near end-side of the Galactic Bar (Clark et al 
2009; Dayiesetal. 2007; Figeretal. 2006; Negueruela et a" 
120111 l2010h . We restricted' our analysis to a sample of 
119 spectroscopically confirmed RSGs in clusters, because 
their association with a c l uster confirms the luminosity class 
(Bernabei & Polcaro 2001; Caron et al. 2003; Clark et al. '20091 
Davies et al. 2009a, 2008; Eggenbergeret al. 2002; FigeretaL 
2006, 1999; Mengel & Tacconi-Garman 2007; Mermifliod et aL 
2008; Messineo et al. 2011, 20 10: .Negueruela et aL>201 Li20ia 
Pierce et al. 2000; Skiff 20l(%. When comparing the distribution 
in colors of the 119 RSG stars in clusters with all RSGs listed by 
Skiff et al., however, no significant differences are seen. 

IVerhoelst et alJ (l2009l) studied the dust properties of 0-rich 
RSGs envelopes and found a dust condensation sequence re- 
sembling that of AG Bs ( Kemper et al. 2001; Speck etal. 200^ 
ISvlvester et anil999l) . Stars with low mass-loss rates are rich of 
oxides and alumina, while high ma ss-loss stars show olivine, and 
crystalline silicates (IWatersll20ip|). Metallic Fe or amorphous 
carbon could also be present (Ver hoelst et alJ 120091: IWatersI 
2010). Buchanan et al. (2009) report PAH emission at 6.2 fim in 
some LMC RSGs. PAHs are not detected in mid-in frared spec- 
tra of Galactic RSG stars obtained with ISO-SWS ( Sloan et all 
2003). 

Among our sample of 119 RSGs, 69 matches were found 
with GLIMPSE counterparts and 79 with MSX sources. 34 stars 
were detected by both the MSX and the GLIMPSE surveys. The 
mean value of their [8.0] - [A] color is +0. 1 1 mag with a standard 
deviation of 0.20 mag. 



2.4. Mass-losing AGB stars 

AGBs are stars with masses below 8 Mq, with a degenerate core 
consisting of C and O, and two burning shells. The inner shell is 
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burning He, while the external one is burning H. The presence of 
two different sources of nuclear energy determines an interplay 
between these two shells, causing thermal pulses with inward 
and outward motion of matter. Carbon can be brought onto the 
stellar surface, and the AGB star is called a C-star, when the 
ratio of C over O is above 1. The relative numbers of O and 
C-rich AGBs strongly depend on metallicity. We only consider 
here O-rich AGBs, since AG Bs in the inner Galaxy are mostly 
O-rich types (ISevensteiill999h . 

Several different names are used to indicate specific sub- 
classes of AGBs, based on their pulsation properties and/or cir- 
cumstellar properties (pulsato rs, masing stars). A detailed review 
on AGB stars is presented bv lHabing & OlofssonI (l2003h . Miras 
are classically defined as AGBs with visual amplitudes larger 
than 2.5 mag; their periods typically range from 150 to 1500 
days. SRs are regular pulsators with visual amplitudes smaller 
that 2.5 mag; their periods range from 35 to 250 days. AGBs 
may have circumstellar envelopes, where maser emission origi- 
nates. AGBs may be called SiO masing stars, or OH/IR stars, if 
they have SiO or OH maser emission. 

Mass-loss in AGBs is caused by stellar pulsation, and in- 
creases going from SR stars to Miras and to OH/IR stars, i.e. 
with increasing periods. It scales with stellar luminosity (M oc 
p-'' in SRs) and may range from ~ 10'^ to ~ 10'-^ Mp yry' (e.: 
lAlardetal.ll200U iHabing & OlofssonI l2003t lOrtiz et alJl2()02 
SR stars have typical mass-loss rates from ~ 10"^ to ~ 10"* M© 
yr"', while Miras have mass-loss rates from ~ 10"^ to ~ lO"'* 
Mq yr"'. A median mass-loss rate of 3 x 10"^ yr"' was esti- 
mated for a sample of OH/ IR stars in the Galactic center region 
(IHabing & Olofssonll2003L and references therein). 

In order to analyze the full range of near- and mid- 
infrared colors of AGBs, we collected samples of AGBs with 
known pulsational properties. About 300 S R stars were de- 
tected in the MACHO and ISOGAL sur veys (lAlard et al.ll200U 
lAlcock etalJI 19991; ISchuller et al1l2003h . We found 66 matches 
with GLIMPSE point sources and 17 with MSX point sources. 
The 9 sources with available ([8.0] - [A]) colors yield a mean 
color of 0.45 mag with a standard deviation of 0.33 mag. 

A sample of 409 large amplitude variables (LAV) was de- 
tected in the centr al 24 x 24 arcmin^ of the Galaxy in TiT-band by 
iGlass et aL I (l200lh . LAVs coiTespond to optical Miras, with pe- 
riods from 150d to 800d. We found 291 counterparts to the 409 
LAVs in the GLIMPSE catalog and 95 in the MSX catalog. For 
the 7 1 LAVs with measurements at ~ 8 /^m from both MSX and 
GLIMPSE surveys, the mean [8.0] - [A] value is 0.24 mag with 
a standard deviation of 0.8 mag. SiO maser was detected in 77 
of these LAVs (54 with a GLIMPSE counterpart and 40 with an 
MSX counterpart, Imai et al. 2002). 

The sample of SiO masing stars of iMessineo et al.l (l2002l) 
consists of 271 stars, which were color selected from the DENIS, 
ISOGAL, 2MASS and MSX catalogs. Information on variabil- 
ity and photometric properties suggests that the whole sam- 
ple consists of Mira-like stars (Messineo et al. 2004). This is 
a sample of mostly AGB stars, but may Ukely includes a 
few RSGs ([M essine 0II2OO4I). Positions have an accuracy of 1" 
(IMessineo et al. 2004). For the SiO masing stars, we found 177 
GLIMPSE and 258 MSX counterparts. The 125 sources with 
both 8jum measurements have a mean ([8.0] - [A]) =0.34 mag 
and a standard deviation of 0.42 mag. 

Along the text we will con sider as an SiO maser sample 
the SiO masing stars by M e ssineo et alJ (l2002l) plus the 77 
masing LAVs by iGlass et ali (1200 ll) ; llmai et all (|2002[) . while 
we will refer to the remaining LAVs bv IGlass et all (I2OOII) 



simply as LAVs. Ma ser and color properties are discussed in 
IVerhevenetd] (120111) . 

OH/IR stars are large amplitude AGBs with periods gener- 
ally longer than 600d (even exceeding 1500 days), and denser 
circumstellar e nvelopes than Mir as or SRs. The sample of 766 
OH/IR stars by ISevenste^ (|2002|) consists of two sets of maser 
stars detected with the VLA and ATCA interferometers (|Z| < 
45° and \h\ < 3.5°). Positional uncertainties are typically 1". 
We found 295 GLIMPSE counterparts and 610 MSX counter- 
parts. 90 OH/IR stars were detected at 8yum by both MSX and 
GLIMPSE. The mean ([8.0] - [A]) is 0.35 mag with a standard 
deviation of 0.52 mag. 

The different numbers of matches found in the GLIMPSE 
and MSX catalogs are due to difference in their sensitivity and 
sky coverage. While MSX covers the whole Galactic plane to 
+5° of latitude, the GLIMPSE survey is confined to a narrow 
latitude range ±1°. OH/IR stars are bright at 8/vm; many were not 
included in the GLIMPSE catalog, because of saturation. The 
MSX detection limit in A-band is 100 mJy ([A] - 6.9 mag), 
while the Spitzer/IRAC 8/im channel saturates at 444 mJy ([A] = 
5.3 mag), and has a 5 sigma detection limit at 2mJy ([8] = 11.2 
mag). 

3. Interstellar extinction 

We adopt a near-infrared extinction pow er law ( Ak ^ 



(Fritz et al. 201 1; Indebetouw et al.l 2005; Messineo et al 



Nishivama et al.n2006; ,Riek e et al."19 85u .Stead & Hoar&.2009. . 
with an index of a = -1.9 (Messine o et al. I I2005I) . This index 
provides consistent measurements of extinction in the versus 
(J-Ks) and Ks versus (H-K^) diagram, also in fields with very 
high extinction Ak,> 2 mag. Furthermore, this index is consis- 
tent with that measured with hydrogen recombination lines by 
iLandiniet all (11984 . A recent historical review of index mea- 
surements is presented by Fritz et alj ( 1201 Ih . 

For the GLIMPSE measurement s, we used the color ex - 
cess ratios E^-ks/Ej-ks calculated bv llndebetouw et al] (l2005l) . 
and we re-calcu la ted th e extinction ratios using Eq n. 1 in 
Indebeto uw et al.l (l2005h . and Aj/A/cs from Messine o et al.l 
(2005). The assumed and resulting quantities are listed i n Table 
[1] They are lower than those first reported by Indebetouw et al.l 
. Lower values are also measured toward the Galactic cen- 
ter dR-itz et alJl20ni) . 

The use of slightly different extinction ratios does not affect 
our analysis, because the average excess colors are much larger 
than the uncertainties due to interstellar extinction. 

4. Color-color diagrams 

In order to identify photometric criteria that might separate 
early-type stars from late-type stars and identify specific types 
of stars (e.g. WRs and RSGs), we analyzed the distribution of 
evolved stars in several color-color diagrams. 



4. 1 . H-Kf, versus J - H diagram 

Because of their different temperatures, "normal" late-type stars 
and early-type stars form two parallel sequences with increasing 
interstellar extinction in the H-K^ versus J - H diagram. This 
diagnostic tool is suited to only classify the bulk of a stellar pop- 
ulation; it relies on the assumption that mass-losing objects are 
rare. Mass-loss changes the stellar energy distribution. 

Dust grains in the circumstellar envelopes of late-type stars 
absorb the stellar light and re-emit at longer wavelengths 
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12 3 4 

H-Ks 



Fig. 1. 2MASS J -H vs . H-K^ diag r am o f Galactic evolved stars. This diagram does not allow for a photorri etric classification. SR 
variables are taken from'Alard et al.' ('200 D, large amplitude varia bles (LAVs ) stars from iGlass et al.l (1200 lb . and AGBs with SiO 



maser emission from Messineo et al. (2002j). OH/IR stars are from Sevensterl (l2002h . The s election of RSGs i s described into the 



text. WR stars are from Messineo et al. (2011, 2009); van der Huchf ('2OO1'), LBVs are from 'C lark et al.1 (l2005h : lGv"aramadze etaP 
( 12010); Mauerhan et al. (2010, 201 1); Messineo et al. (2011J. The continuous line near the bottom left corner is the locus of dwarf 
stars ( Koornneel 1983h . two dashed Unes mark the traces of an 09 (lower) and an M5 (upper) star with increasing interstellar 



extinction. A reddening vector for Ak, = 1 mag is also shown dMessineo et alJl2005l) . Typical errors are within 0.05 mag in both axis 

Table 1. Infrared extinction ratios (A[^]/Ak,) for the 2MASS and SPITZER/IRAC filters. 

AjIAks AhIAk, Ap.6]/A;f, A[4.5]/Akj A[5.8]/As-j A[8.o]/Aki Reference 

2.86 1.66 Messineo et al. (2005) 

2.5 1.56 Rieke et al. (1985) 

2.5+0.2 1.55 ±0.1 0.56 ±0.06 0.43 ± 0.08 0.43 ±0.10 0.43 ±0.10 Indebetouw et al. (2005) 

2.86 1.66 0.44 0.31 0.29 0.31 Indebetouw et al. (2005) and Messineo et al. (2005 



(e.g. iGroene we gen &de.TongllI993h . In the solar neighborhood, 
AGBs with low mass-loss rates < 10"^ Mq yr"' have de- 
reddened J-K^ colors, {J -K)o, between 1.2 and 1.6 mag, while 
AGB wit h mass-loss frorn IQ-^ - 1Q-* have {J - K)o from 2 to 
6.5 mag dMessineo et al ] |2005L and references therein). AGBs 
can have up to 40 mag of visual extinction. Figure [1] shows the 
H-Ks versus diagram of our samples of evolved stars. Two 
continuous lines are plotted, showing the loci of a "naked" M5 
star (left/higher d ashed line) and of an 09 star (right/lower line) 
(lKoornneedl983h wit h varying in terstellar extinction, as well as 
the reddening vector dMessineo e t al. 2005). It is not possible to 
classify stars with this diagram alone, e.g. all Miras would be se- 
lected as OB stars. RSGs, SRs, and Galactic center LAV stars ap- 
pear located between these two lines. SiO and OH/IR stars show 
even redder H-K^ colors, which are due to continuum absorp- 
tion in //-band by water (e.g.. Fig. 6 and 7 in fe lum et al. 2003; 
Comeron et aP 12004 iF^gel & Whitfordlll987i: IMessineo et alJ 
20051k 



WRs seem to define a sequence in Fig. [T] which lies be- 
low the sequence of "normal" OB stars (H-K^ from 0.0 to 1.5 
mag and J-K^ from 0.0 to 2.0 mag). The distribution of WRs 
overlaps that of evolved late-type stars (OH/IR and SiO mas- 
ing stars). Different physical mechanisms can generate similar 
H-Ks and J - H colors; LBVs have J—K^ colors from ~ 1.0 
to ~ 5.0 mag, and H-K^ from ~ 0.2 mag to ~ 2.0 mag. 
A sample of five LBVs (AFGL2298, Pistol Star, WR102ka, 
LBV1806-20, GO.120-0.048, FMM362) have redder J-K,, 
which are consistent with distant sources at large interstellar 
extin ction (Ak,=3.15 2.99 2.78, 3.0, 3.26, 3.42 mag, respec- 
tively'Barniske et al."2008|; iBibbv et al.ll20oi IClark et aUlIool 
Mauerhan et al...2010.) . 



Since circumstellar absorption moves the stars on the di- 
agram almost along the direction of the interstellar reddening 
vector, it is impossible to distinguish between circumstellar and 
interstellar reddening with this diagram. An estimate of total (in- 
terstellar plus circumstellar) extinction can be obtained from the 
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Fig. 2. 2MASS/GLIMPSE /rs-[8.0] v ersus J-K^ diagram. S ymb ols are as in Fig. [11 The arrow represents the reddening vector 
following the extinction ratios given bv lMessineo et al.l (l2005h and llndebetouw et al.l (l2005l) . Typical errors are within 0.07 mag in 
both axis. 
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Fig. 3. GLIMPSE [3.6] - [8.0] versus [3.6] - [4.5] diagra m. Symbols are as in Fig . \T\ The arrow indicates the direction of the 
reddening vector following the extinction ratios derived by llndebetouw et all (l2005h . Typical errors are within 0.08 mag in both 



axis. 
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reddening in J-K^, Akv ats = E(( J-Ks)) * 0.54 (iMessineo et alJ Table 2. Ql values for normal stars (without mass-loss). 
I2005h . 



4.2. ^rs-[8.0] versus J-K^ diagram 

A useful diagnostic for mass-loss is a K^-X versus J-Kf^ 
diagram, where X is a mid-infrared magnitude, e.g. the 15 
fim measurements made by ISOGAL, o r the 8 ;um measure- 
ments made by MSX and GLIMPSE (e.g. IMessineo etal.ll2005t 
Ivan Loon et al E00l.InFig.l2l we show a/rs-[8.0] versus J-K^ 
diagram. The reddening vector is determined u sing the near- 
infrared e xtinction law by | Messineo et al.l (l2005l) . and the color 
excess bv llndebetouw et aL (12005 ) (see Sect. 4). A separation 
between interstellar and circumstellar extinction is detected in 
this diagram, since increasing circumstellar reddening do not 
move the star along the direction of the interstellar reddening 
vector Deviations from the reddening vector provide evidence 
for mass-loss (e.g. IMessineo et al.l2005l , and references therein). 

SRs are almost located on the reddening vector, while OH/IR 
stars show large /rs-[8.0] color excesses (up to 8.0 mag). There 
is a large overlap between the colors of AGBs (OH/IR stars, SiO 
masing stars, and Miras) and RSGs stars. 

WRs follow a separate seq uence in the A's-[8.0] versus J-Kf^ 
diagram (Hadfiel d"etal]|2007l) . This WR sequence suffers con- 
tamination from late-type stars. 



4.3. [3.6] - [8.0] versus [3.6] - [4.5] diagram 

In Fig. [3] we show a [3.6] - [8.0] versus [3. 6] - [4.5] diagram 
of evo lved stars. The reddening vector from Indebetouw et alJ 
(l200l is also marked. Evolved stars do not follow the interstel- 
lar reddening vector, but, for a given [3.6] - [4.5] color, they dis- 
play redder [3.6] - [8.0] colors. We expect interstellar extinction 
in the IRAC bands below 0.3 9 mag even in the central regions of 
the Milky Way (Ak~ 3 mag) (llndebetouw e t al. 2005). Since the 
ranges of colors seen in Fig. |3] are much larger than this value, 
we can extract indication on intrinsic colors for each stellar type. 

In Fig. [3] AGBs (SRs, Miras, and masing stars) form a se- 
quence of increasing [3.6] - [8.0] colors, from 0.0 (SRs) to 5.0 
mag (OH/IR stars). This is likely due to a combination of in- 
creasing mass-loss rates and water absorption. SR and LAVs 
stars display increasing [3.6] - [4.5] colors, with increasing 
[3.6] -[8.0] colors, likely due to increasing strengths of water ab- 
sorption at ~ 3fj.m. SiO masing stars display peculiar [3.6] - [4.5] 
colors, bluer than those of SRs and LAVs (without SiO masers). 
A possible cause is a difference in the water strengths and pos- 
sible SiO absorption at ~ 4 //m and CO2 absorption at ~ 4.3 
/Lim dCamii 120021) . The distribution of RSGs resembles that of 
masing AGE stars in this plane. Some RSGs have the bluest 
[3.6] - [4.5] colors than any other stellar type, as alrea dy noted 
in the LMC, likely due to CO (e.g.'Yang & Jiang"201l'). 

Few ISO-SWS spectra of late-type'stars (Sloan et al1l 2003h 
are shown in Appendix. 

The bulk of WR stars have [3.6] - [8.0] from 0.5 to 1.5 mag 
and [3.6] - [4.5] from 0.0 to 0.7 mag. The colors of early type 
stars, WRs and LBVs, overlap with those of late type stars. 

5. The Q\ and Q2 parameters 

5.1. The Ql parameter 

The Ql parameter, Ql = (J - H) - 1.8 x (H-K,), is a mea- 
sure of the deviation from the reddening vector in the H-K^ 
versus J-K^ plane, following the infrared extinction law by 



Sp. type 


I 


Ql 


III Ql 


V Ql 




09 


I 


-0.02 


III 


V -O.I 




AO 


I 


0.07 


III 


V 0.01 




FO 


I 


0.13 


III 


V 0.13 




G3 


I 


0.31 


III 0.36 


V 0.23 




KO 


I 


0.36 


III 0.41 


V 0.29 




M5 


I 


0.72 


III 0.68 


V 0.38 




Table 3. 


Typical Ql and Q2 ran 


ges for the considerec 


of stars. 












Type 


Ql (min) 


Ql (max) 


Q2 (min) 


Q2 (max) 


RSG 




-0.06 


0.44 


-3.38 


1. 12 


WR 




-1.19 


O.ll 


-5.02 


-1.02 


LBV 




-0.34 


0.06 


-2.23 


-0.73 


OH/IR 




-1.56 


0.44 


-12.99 


-0.49 


SIO 




-0.60 


0.30 


-7.31 


0.69 


LAV 




-0.17 


0.53 


-3.01 


1.49 


SR 




0.02 


0.52 


-1.04 


1.46 



IMessineo et alj (|2005|) . This quantity does not depend on in- 
terstellar extinction. The Ql parameter was first introduced to 
photometrically se lect counterparts of high-mass X-ray binaries 
dNegueruela & Sch urch 200'^; it has been successfully used for 
searching clusters of RSGs (e.g. Clark et al. 2009; Davies et alj 
2007; Figeretal. 2006; Negueruela et al. 20n, 2010). Since 
stars with different spectral types have different intrinsic colors, 
they will also have a d ifferent value of <9 1, as shown in Table |2] 
with intrinsic colors bv lKoornneef) d 1 9831) . Ql is typically around 
zero for "n ormal" OB stars, while it is around 0.4 mag for K- 
giant stars (iKoornnee f 1983). 

An infrared excess due to free-free emission in the winds of 
hot stars (e.g. WRs), or to dusty circumstellar envelopes (e.g. 
AGBs or RSGs), can modify the Ql parameter of stars. Mass- 
losing late-type stars may have negative Ql values (similar to 
those of early type stars). The Ql parameter alone does not allow 
distinguishing between early and late-type stars. 

In Fig. |4] the Ql parameter is plotted against the A's-[8.0] 
color. Two panels are shown to separately plot the Ql versus the 
/:s-[8.0] colors from the 2MASS and GLIMPSE catalogs, and 
the Ql versus K^-A colors from the 2MASS and MSX catalogs. 
The numbers of stars per type are different in the two panels. 
Most of the OH/IR stars are detected by MSX, but are saturated 
in GLIMPSE. Galactic center SRs could only be detected with 
the more sensitive GLIMPSE survey. Typically, Ql is less than 
0.1 mag for windy early-type stars, and ranges from -0.1 to 0.4 
mag for the bulk of known RSG stars. AGB SRs and AGB LAVs 
have similar ranges of Ql values, which are also similar to those 
of RSGs. AGBs with SiO or OH maser emission have broader 
Ql ranges, OH/IR stars have values down to ~ -1.6 mag (see 
TableO. Late-type stars move to redder /Tj- [8.0] with decreas- 
ing Ql values. Most of the WRs have negative Ql values. LBVs 
have Ql values range from 0.1 to -0.35 mag. 

In Fig. |4]the Ql values are interstellar extinction indepen- 
dent, while the /rs-[8.0] and K^-A values are affected by inter- 
stellar reddening (E(/rs-[8.0])), and, therefore, interstellar red- 
dening can only move a data point on a horizontal line. We 
detected a diagonal cut-off in the distribution of Ql versus 
/rs-[8.0] values of late-type stars in Fig.|4] This is unrelated to 
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Fig. 4. Ql parameter versus A's-[8.0]. Symbols are as in Fig.[Tl In the top panel the GLIMPSE 8 //m measurements are plotted, 
while in the bottom panel the MSX 8 um measurem ents (A-band) are plotted. Th e arrow represents the reddening vector following 
the extinction ratios given by'Messineo et al.' ('2005') and Indebetouw et al. (2005). The dashed lines are a linear fit to the Ql versus 
the dereddened (/rs-[8.0]) values of SiO masing stars (,MessineQy2004i) for Ak=0,1,2 mag (from bottom to top). Typical errors are 
0.05 mag in the x-axis and 0. 14 mag in the y-axis. 



interstellar extinction, and arises from the intrinsic stellar energy tion significantly changes with increasing mass-loss also in the 
distribution (SED) of late type stains. The stellar energy distribu- J, H, and Ks bands. Bolometric connections in /Tj-band (BCks) 
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and in A-band (BCa) as a function of Ks-A were calculated for 
a sample of SiO masing stars in the inner Galaxy, and compared 
with those of AGBs in the solar neighb ours (Chapters, Appendix 
A2 of Messineo 2004; Olivier et al.i l200U IWhitelock & Feasli 
[1994 .Whitelock et aL .2000. ^003). BCk^ decreases with in- 
creasing Ks-A, and for Ks-A=6 mag BCk^^ 0.8 mag, i.e. ~ 2.3 
mag smaller than values for "naked" M-type stars (2.9-3.1 mag) 
(e.g. Blum et al. 2003; Frogel & Whitford 1987). 

We obtained the following relations bet ween the 01 an d 
A's-[8.0] values of SiO masing Mira-like stars dMessinec l2004l) : 

Ql = 0.50(+0.17) - 0.22(+0.05) x BCa 

Ql = 0.34(±0.05) - 0.20(±0.04) x (Ks - A)o. where K^-A 
ranges from ~ to ~ 4 mag. 

These relations prove that changes in the SEDs cause vari- 
ations in the Ql values, and generate the cut-off seen in Fig.H] 
Reddening in K^-A, due to interstellar extinction, and stellar pul- 
sation smooth such a cut-off by smearing it over K^-A =~ 1.5 
mag. For an Ak, of 1, 2, and 3 mag, E(/:s-[8.0]) =~ 0.31, 
~ 0.62, and ~ 0.93 mag, respectively. This relation nicely ex- 
plains the diagonal cut-off observed in Fig.lD OH/IR stars with 
much redder colors (K^-A > 4.0 mag) may fall below the rela- 
tion found for SiO masing stars. 

This relation between Ql values and A's-[8.0] colors allows 
a separation of mass-losing late-type stars from hot windy stars 
(e.g. WR stars). A full revision of bolometric corrections in sev- 
eral near and mid-infrared bands will be the topic of a forthcom- 
ing paper. 

A few peculiar WRs fall above these Unes; WR 125 has a 
Ql value 0.8 mag above this line when using its MSX 8 fim 
measurement. This so urce is known to be dusty and variable 
(I Williams et al.l |1994|) . LBVs cannot be identified in this dia- 
gram. Their (Ks -A)o colors expand over the region of WR stars, 
as well as over that of late-type stars. 

5.2. The Q2 parameter 

We define a new parameter Q2 = {J-K^) - 2.69x(Ks-[8.Q]), 
which is a measure of the deviation from the reddening vector in 
the (7-^s) versus .^s-TS-OI plane , following the infrared extinc- 
tion law of'Messineo et al.' ^2005) and mid-infrared color excess 
ratios by Indebetouw et al. (2005). Q2 can be thought of as a 
measure of an excess that is only due to a circumstellar shell 
(E(A's-[8.0])j/,e//), and is independent of interstellar extinction: 

Q2^(E(Ks-[m)sheid- 

We prefer to use the excess in the (J-Ks) versus (A's-[8.0]), 
rather than in the {H-Ks) versus (A's-[8.0]) plane, because 
AGBs may have strong absorption in H-band by gaseous wa- 
ter (see Fig, ffl or Blumetal. 2003 ; .Frogel & Whitford 1987; 
iMessineo etalllloash . 

The Q2 parameter can be used to statistically select stars 
with infrared excess. Such excess can be due to dust emission 
(late type stars), or free-free emission in shocked winds (WRs); 
for LBVs it is generally due to free-free emission plus dust emis- 
sion. 

In Fig.|5] the Q2 parameter is plotted against the [3.6] - [4.5] 
color. Q2 is a parameter independent of interstellar extinction, 
interstellar reddening in the [3.6] - [4.5] color are small because 
the e xtinction ratios A[3.6i/Ak , =0.44 and A[4 5]/Ak, =0.31 are sim- 
ilar jlndebetouw et al.ll2005h . The reddening in [3.6] - [4.5] is 
0.13, 0.26, and 0.39 mag, for Aks = 1, 2, 3 mag, respectively. 

The samples of masing stars (which are mostly Mira-like 

stars) show typical Q2 values from 13 mag to ~ 0.7 mag, but 

mostly smaller than -1.0 mag. SR stars and RSGs have mostly 



Q2 between 1.5 and -1.1 mag. By combining color and magni- 
tude information it is possible to select bonafide candidate RSG 
stars, since about 70% of the stars in the range of Q2 from -1.1 
and 1.5 mag, and [3.6] - [4.5] < -0.4 mag, are RSGs. WR stars 
have Q2 values from -5.0 mag to -1.0. LBVs have a narrow 
range of Q2 values from -2.2 to -0.7 mag. 

5.3. Ql versus Q2 

In Fig. |6] we plot the Ql versus Q2 values for the samples of 
known evolved stars described in Sect. 3. We estimated the Ql 
and Q2 ranges of Table [3] by making histograms of their values 
and retaining only bins with at least 4 elements; the enclosed 
fractions are generally above 90%. 

Early type stars (with the exception of LBVs) and late type 
stars are distributed differently in the Ql versus Q2 diagram. 
Early-type stars are mostly located above the relation Q2 = 
1 1.25* gl -2.38 mag, while mass-losing late-types (with the ex- 
ception of SRs) are mostly located below it. However, a clearer 
separation between early and late-type stars appears in the Ql 
vs K - A diagram (see Fig.|4|, where early and late-type mass- 
losing stars can be separated with a 20% error (see Table |4|i. 
Specific searches can be optimized by adopting more stringent 
constraints on Ql and Q2, and by using luminosity information. 

AGBs form a sequence of decreasing Ql and Q2 values go- 
ing from SRs to OH/IR stars. The bulk of SRs have Ql from 
0.0 to -1-0.5 mag, and Q2 from -1.0 to 1.5 mag. LAV stars have 
similar Ql values, but their Q2 values have a larger range (from 
-3.0 to 1.5 mag). Masing stars have a slightly larger range of Ql, 
from -1.6 mag to 0.4 mag, and generally negative Q2, down to 
-13.0 mag. 

The Ql and Q2 values of RSGs are similar to those of LAV 
stars. Their Ql values range from -0.0 to 0.4 mag, and Q2 are 
mostly from -3.4 to 1.1 mag. A large fraction (56 %) of the 
RSGs appear concentrated in the region 0.1 < Ql < 0.5 mag 
and -1.1 < Q2 < 1.5 mag, and 26 % in the narrower region 
0.1 < Ql < 0.5 mag and 0.5 < 22 < 1.5 mag. The fraction 
of LAVs and masing stars falling in the latter region is below 
2%. SR stars cannot be distinguished from RSGs on the basis of 
colors only (49% of SRs fall into the narrow region). Since SRs 
are intrinsically much fainter, the degeneracy can be eliminated 
with luminosity information. 

WR stars have distinct Ql and Q2 values; typical Ql val- 
ues are smaller than ~ 0.1 mag and Q2 values smaller than -1.0 
mag. LBVs have typical Ql values from -0.3 to 0.0 mag, and 
Q2 values from -2.2 to -0.7 mag. A distinct group of three 
LBVs (Wra751/IRAS1 1065-6026, AFGL2298, and HD168625) 
appears at Ql ~ -0.2 mag and Q2 ~ -9.0 mag. HD 
168625 and AFGL2298 have large Q2 values because of their 
strong PAH emission lines (Pasauali et al.l l2002l iPeeters et al.l 
2002; Umanaetal. 2010). This suggests that the spectrum of 
Wra75 1/IRAS 1 1 065-6026 is also dominated by strong PAHs. 



6. Suggested selection criteria 

There is overlap among the various samples in the color-color 
diagrams shown in Figs.[ri|2] and|3] It is not possible to uniquely 
define a photometric classification scheme that is only based on 
these diagrams. Furthermore, each color depends on interstellar 
extinction. 

We propose a new classification based on two extinction free 
parameters, Ql and Q2 (see Figs.|4]|5] and|6]l, in which we iden- 
tify color-windows mostly populated by RSG stars and massive 
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Fig. 5. Q2 parameter versus the GLIM PSE [3.61 - [4.51 color. S ymbols are as in Fig. [T] The arrow represents the reddening vector 
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windy hot stars (e.g. WRs). Contamination by other stars can be 
reduced by taking into account luminosity and distance informa- 
tion. 

The Q2 values is based on near and mid-infrared data, which 
are not taken simultaneously. In order to reduce the uncertainty 
due to non simultaneity, we use information on variability, e.g., 
flags provided in the MSX catalog, comparis on of the 8 um 
measurements from GLIMPSE and MSX (e.g. iRobitaiUe et alj 
[2007), and of the J and K g measurements fr om 2MAS S and 
DENIS Ks (e.g. lMessineoet al. 2004; Schultheis & Glass 2001). 

In the following, we list our new proposed classification 
steps. 

- With available J,H,Ks magnitudes from 2MASS and 8.0yum 
magnitudes from GLIMPSE or MSX, we calculate the Ql 
and Q2 parameters. 

- We classify as late-types those stars located above the rela- 
tion gl = -0.20x(/:s-[8.0])-H0.34mag(seeFig.|4|,andwe 
classify as early-type stars those stars below it. This relation 
applies only to stars with A's-[8.0] < 4 mag. 

- To explore the impact of each constraint, we select can- 
didate early-type stars with free-free emission in several 
manners. We call "free-1" stars those early-type stars with 
Q2 < -I mag, and "free-2" stars with the additional con- 
dition Q2 > 11.25 xQl - 2.38 mag. We cafl "free-3" stars 
with ext ra additional condition s , the GLIMPSE criteria spec- 
ified by iHadfield et af] (l2007h : iMauerhan et al] (1201 lb. We 
compare the new criteria with those used by Mauerhan et alJ 
(hereafter, called "free-MVMl 1" 201 1). 

- We select late-type stars with 0.1 < Ql < 0.5 mag and 
-1.1 < g2 < 1.5 as candidate RSGs (cRSGl), or more strin- 



gently with 0.1 < gl < 0.5 mag and +0.5 < Q2 < 1.5 mag 
(cRSG2). When using the additional condition [3.6] - [4.5] < 
-0.4 mag, a detection rate of 70% is expected (see Fig. |5]l. 
We call this more restricted color selection as cRSG3. We do 
not consider candidate RSGs with an indication of variabil- 
ity. 

- The degeneracy between candidate RSGs and AGB-SRs can 
be removed by additional information on luminosity, and 
therefore extinction and distance. We suggest to calculate the 
to tal extinction (c i rcums tellar and interstellar) as described 
in iMessineo eTa and to use the model of Galactic 

dust distribution by Drimmel et alj ( l2003h for estimating dis- 
tances. With bolometric corrections, distances, and extinc- 
tion estimates, luminosities can be derived. 

For each of the samples, we ignored the a priori knowledge 
on spectral types, and calculated the fractions of retrieved free- 
free emitters and RSGs, by using the categories defined above. 
The results of our experiments are listed in Table|4] 

The selecti on of early and late- type stars retrieves 85% of 
WR stars from Ivan derHuchtl (1200 lb but only 63% of S R stars, 
and 48% of the OH/IR stars listed bv lAlardetal] (l200lh . 

The Q2 parameter is a good discriminant for selecting free- 
free emitters ("free-1"); additional constraints on GLIMPSE col- 
ors are needed to reduce contamination by late-type stars ("free- 
3"). The pre-selection of early and late-type stars allows retrieval 
of 75% of free-free emitters, but strongly reduces the contamina- 
tion by AGBs; the number of AGB stars erroneously retrieved as 
free-free emitters decreases from 33% ("free-l") to a few percent 
("free-3"). For a comparison with Mauerhan et al. (201 j), w e 
also use their selection criteria (Fig. 1 in .Mauerhan et aljl201 ih . 
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Table 4. Fractions of retrieved types. 



Type 


Early-type 


free- 1 


free-2 


free-3 


free-MVMll 


Late-type 


cRSGl 


cRSG2 


cRSG3 


WRall 


0.86 


0.85 


0.85 


0.78 


0.80 


0.14 


0.02 


0.00 


0.00 


WR MVMl 


0.75 


0.73 


0.73 


0.72 


0.73 


0.25 


0.00 


0.00 


0.00 


WR SFZl 


0.90 


0.90 


0.90 


0.90 


0.90 


0.10 


0.05 


0.00 


0.00 


LBV 


0.67 


0.67 


0.67 


0.33 


0.44 


0.33 


0.11 


0.00 


0.00 


RSG 


0.11 


0.00 


0.00 


0.00 


0.00 


0.89 


0.72 


0.42 


0.22 


OH/IR 


0.52 


0.33 


0.15 


0.00 


0.15 


0.48 


0.04 


0.00 


0.00 


SiO 


0.25 


0.16 


0.16 


0.01 


0.19 


0.75 


0.04 


0.00 


0.00 


AGB-LAV 


0.03 


0.01 


0.01 


0.01 


0.03 


0.97 


0.36 


0.06 


0.00 


AGB-SR 


0.37 


0.02 


0.02 


0.00 


0.00 


0.63 


0.58 


0.35 


0.00 



Notes. For each sample of star, we list the fraction of retrieved early-type (Early-type), three estimates of the fractions of free-free emitters (free-1, 
free-2, free-3, and free-MVMll), the fraction of retrieved late-type stars (Late-type), and three estimates of RSG stars (cRSGl, cRSG2, cRSG3, 
and cRSG4). Only stars with counterparts in the three 2MASS and four GLI MPSE ban ds are consid ered. F or WR stars, we consider the whol e 
sample alltogether (WR all), as well as the newly detected samples by Mauer han et alj (WR MVMl 1 1201 ih and lShara etaP (WR SFZl 1 1201 ih . 
The selection of free-free emitters is applied only to early-type stars; the selection of RSG stars is applied only to late-type stars. 



which we label as free-MVMl 1. The Q1-Q2 based criteria pro- 
vide a detection efficiency comparable to that by Mauerhan et al. 
on the whole sample of WRs. These aut hors applied additiona l 
criteria than the color-color criteria by iHadfield et alJ (l2007h . 
based on the Ks versus J-Ks diagram. We also reclassified sepa- 
rately their sample of 61 WRs, which is composed of more dis- 
tant and obscured WR stars with our classification scheme, and 
retrieved 75% of them as early-type stars, and 73% as free-free 
emitters (free-1). 

The selection of a narrow window in the Ql versus Q2 plane 
with the additional limits on GLIMPSE colors allows selecting 
bonafide RSG stars. 



7. Summary and discussion 

7.1. Summary 

We analyzed 2MASS and GLIMPSE color properties of samples 
of Galactic evolved stars. Samples of O-rich AGB stars (SRs, 
Miras, OH/IR stars, and SiO masing stars), RSG stars, WRs, 
and LBVs were collected from existing literature. Several color- 
color diagrams were analyzed, aiming to identify the best com- 
bination for stellar photometric classifications. 

The J - H versus H-Ks diagram, and the J-K^ versus 
A's-[8.0] diagram are particular useful. In these planes, the col- 
ors of WR stars and mass-losing late type stars deviate from the 
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interstellar reddening vector. WRs have a mid-infrared excess 
due to their free-free emission. Late-type stars often have mid- 
infrared excess due to their dusty envelopes. The deviation from 
the reddening vector in the A's-[8.0] versus J-K^ plane increases 
going from SRs to Miras and OH/IR stars, i.e. with increasing 
mass-loss rates. 

The J - H versus H-K^, and the J-K^ versus /rs-[8.0] dia- 
grams together can be used to discriminate between interstellar 
and circumstellar reddening, and identify different types of stars 
with envelopes. Stellar classification is hampered by color over- 
laps between sources of different types, and scatter due to vari- 
ability and non-simultaneity of the observations. Special color 
windows have been identified for locating WR and RSG stars. 

We investigated two extinction free parameters, Ql and Q2. 
gl is a measure of the deviation from the reddening vector in the 
J -H versus H-K^ plane (Negueruela & Schurch 2007). Q2 is a 
newly defined parameter, which measures the deviation from the 
reddening vector in the J-K^ versus A's-[8.0] plane. The Q\ and 
Q2 parameters allow for an efficient selection of stars with free- 
free emission (e.g. WR stars), and candidate RSG stars. Stars 
with Q\ < -1.0 mag and g2 > 11.25 x Q\ -2.38 mag ai-e can- 
didate hot massive stars with free-free radiation. A large number 
(~ 40%) of Galactic RSGs are found with 0.1 < gl < 0.5 mag, 
and -1.1 < Q2 < 1.5 mag. Selections can be further improved 
with additional conditions (e.g. with the [3.6] - [4.5] color). 

A combination of mid- and near-infrared measurements al- 
lows to statistically distinguish between evolved mass-losing 
stars and early-type stars. 

7.2. Selection of free-free emitters 

We have found a clear separation between massive stars with 
free-free emission (e.g. WR stars) and late-type stars in the Ql 
versus ^rs-[8.0] diagram, and in the Ql versus Q2 diagram. 

A combination of [3.6] - [8.0] versus [3.6] - [4.5] and 
A"s -r8.0] versus J-K^ diagrams (Figs. |3] and |2]l were chosen 
by iHadfield et"an (12007.) to best iden t ify W R stai-s from other 
field stars. Recently, Mauerhan et al.l (1201 Ih revised the selec- 
tion of WR stars with additional constrains on the H-K^ ver- 
sus J-K^ diagram, and versus J-K^ diagram, to obtained a 
spectroscopic detection rate of 95% in ear ly-type stars, a n d 20% 
in WR stars. The criteria suggested by IHadfield et al.l (l2007h 
and lMauerhan et al] ( 1201 Ih depend on interstellar extinction, and 
suffer some contamination by dust enshrouded AGBs. Our new 
selection based on 21 and Q2 has the advantage of simplicity, 
based on a single diagram that is independent of interstellar ex- 
tinction. A fraction of WRs (about 15 % ), however, are found 
in the region of late-type stars, and are missed by our selection 
of early type stars - among them several WC stars, which are 
known to have dusty envelopes. 

The complex structures of LBV envelopes, which are dusty 
and often extended, result in a broad range of colors. Q2 varies 
from -2.2 to -0.7 mag and Ql from -0.34 to 0.0 mag. LBVs 
occupy a broad color spaces, and cannot be identified with pho- 
tometry of point sources alone. 



7.3. Pulsation and stellar colors 

A correlation is newly found betwee n pulsation types and 
GLIMPSE colors. lAlvarez et al.l (12000') writes: "Pulsation pro- 
duces much more extended atmospheres, and in addition dense 
cool layers may result from the periodic outwards running 
shocks. In various ways pulsation thus leads to the existence 



of regions where relatively low temperatures are combined with 
relatively high densities, conditions that favor the formation of 
water, and dust." 

For AGBs, we have located a sequence of increasing [3.6] - 
[8.0] colors with increasing pulsation amplitudes (from SRs to 
OH/IR stars). The behaviors of the [3.6] - [4.5] color with spec- 
tral types is more complex. Generally, SRs and LAVs (without 
SiO maser) show a nice correlation of increasing [3.6] - [4.5] 
colors with increasing [3.6] - [8.0] colors, but the subgroup of 
SiO masing stars appears to be more scattered and to have bluer 
[3.6] - [4.5] colors. In the Appendix, a few ISO-SWS repre- 
sentative spectra of evolved late-type stars are shown. The in- 
crease in [3.6] - [4.5] colors going from SRs to LAVs is due to 
strong continuum absorption by water in the GLIMPSE 3.6 yum 
band (Matsuura et al. 2002). In contrast, the bluer color of mas- 
ing AGB stars is likely due to a combination of stronger water 
absorption (3.6 yum band), and presence of absorption due to SiO 
and CO2 molecules in the 4.5 fim band. 

Several parameters, such as luminosity, rotation, turbulence, 
and metallicity, play an important role in determining the dust 
chemistry of RSG envelopes. RSG stars with regular pulsation 
display redder [3.6] - [4.5] color than irregular RSGs, as sug- 
gested by existing ISO-SWS spectra (see Appendix). Pulsation 
appears, therefore, a key parameter for mass-loss rates also in 
RSG stars. The GLIMPSE [3.6] - [4.5] and [3.6] - [8.0] colors 
of RSGs correlate with mass-loss rates (see Appendix). RSGs 
show a broad range of [3.6] - [4.5] colors, and can be bluer than 
AGB stars. Photometric monitoring of Galactic RSGs, as well as 
follow-up mid-infrared spectroscopic observations of the 3 fim 
region, is needed to further explore a correlation between chem- 
istry and pulsation properties in RSGs, as suggested by Fig. [3] 
The RSGCs share the sa me location in the G alaxy, and likely 
have similar metallicity ( Davies et al.l l2009bh . but their RSG 
members have a large spread in [3.6] - [4.5] colors. This sug- 
gests that pulsation affect the mid-infrared properties of RSGs 
more than metallicity. 



7.4. Selection of RSG stars 

We confirm that the se l ection of RSGs based o n Ql v alues 
adopted by Clark etal. (l2009h . iNegueruela et alJ (120 lOl) . and 
iNegueruela et al.l (l201 1 ) (as remarked already in these works) 
excludes possible dust enshrouded RSG stars, and its strongly 
contaminated by AGBs (SRs and Miras). Selection of bonafide 
RSG stars solely based on the Ql parameter is impossible. The 
Ql parameter remains, however, a valid tool to select clus- 
ters with RSGs. RSGs are often found in massive c lusters, be- 
cause of their young ages (< 30 Myr) (e.g. iDavies e t al. 2002 
Figer et al. 20 061). Masing AGBs are typically found in isolation 
(e^.Messineo et alJ 120021 : ISevensten l2002t ISjouwerman et al.1 
\l9m. 



We propose a selection of candidate RSG stars based on a 
larger set of constraints: the Ql and Q2 values, variability infor- 
mation, [3.6] - [4.5] colors, and apparent magnitudes. We can 
distinguish between interstellar and circumstellar extinction, us- 
ing the Ql and Q2 parameters. First distance estimates can be 
obtained by assuming a n interstellar-extinction versus distance 
relation (Drimmel et al.ll2003h : stellar luminosities can, then, be 
estimated by using adeguate bolometric corrections. 
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7.5. Final remarks 

While our new criteria do not allow classification of all IR bright 
Galactic stars detected by GLIMPSE or MSX, the number of 
known WR stars and RSGs can highly be increased. The pro- 
posed criteria are particularly useful for studying stellar popu- 
lations, for which an average extinction and distance can be as- 
sumed. They will allow to efficiently detect massive stars in giant 
HII regions and supernova remnants. 

The presented datasets open new possibilities for address- 
ing, and pose questions concerning our understanding of stellar 
reddening, bolometric corrections, and extinction ratios. 

Appendix A: ISO-SWS spectra of late-type stars 

Molecular absorption and/or emission and dust features affect 
the [3.6] - [4.5] and [3.6] - [8.0] colors. In order to understand 
the observed trend of increasing [3.6] - [4.5] and [3.6] - [8.0] 
colors of late-type stars with pulsation types, we analysed ISO- 
SWS spectra (see Fig. [B- From the flux-calibrated ISO-SWS 
library of Sloan et al. (2003), we selected spectra of stars with 
known spectral types and variabilities: AGB-Miras, AGB-SRs, 
and OH/IR stars, and RSGs. 

Water, SiO, CO2, and CO absorption may be present in the 
3.0 to 5.0 yum spectra l region of late-typ e stars (Cami' '2002^, 
iMatsuura et al. 2002; Sylvester etalJlT999). ISO-SWS spectra 
of RSGs and AGB-SRs are remarkably different from those of 
AGB Miras (see Fig.Q]). Typically RSGs and AGB-SRs do not 
have strong molecular signatures in between 3 and 4 jum, and dis- 
play a decreasing flux in this region (GLIMPSE 3.6 fj.m). AGB- 
Miras are dominated by water, SiO, CO2, and CO absorption, 
and the resulting spectral flux between 3 and 4 fim can have a 
positive slope. This is attributed to a different "MOLspheres", 
layers of molecular material in an extende d atmosphere (Miras) , 
where water is the main source of opaci ty jMatsuura et al.ll2002t 
iTsuii et all [1991 IVerhoelst et al.ll2009h . The st rength of water 
absor ption varies with stellar pulsation phase (iMatsuura et alJ 
l200l. 

IVerhoelst et alj (l2009h explains the continuum shape of 
RSGs in the 2.5-3 /um regions with an extra dust component 
(metallic Fe, amorphous C, or O-rich micron-sized grains) or 
with chromospheric emission. The similarity of the continuum 
of AGB SRs and RSG stars in the 2.5 and 3 jum region suggests 
that the attribution of extra continuum emission to free-free radi- 
ation from chromospheric activities is less likely. Regularly pul- 
sating RSGs have larger mass-loss (as suggested by the strength 
of the silicate emission at 9.7 fim) than RSGs with irregular pul- 
sations (Fig.rri). 

The ISO-SWS flux densities at the effective wavelength of 
the GLIMPSE filters were measured, and flux densities were 
converted into magnitudes using the GLIMPSE zero-points. The 
synthetic [3.6]-[4.5] versus [3.6]-[8.0] diagram is shown in Fig. 
[72] and well reproduces the sequence of increasing [3.6] - [4.5] 
with [3. 6] -[8.0] colors. A clear sequence of redder colors is seen 
in AGB stars, going from SRs, to Miras and OH/IR stars. An in- 
crease in the [3.6] - [8.0] color is mostly due to an increase of sil- 
icate emission at 9.7 jum, which for optically thin envelopes well 
correlates with mass-loss. Their [3.6] - [4.5] colors are redder 
than those of AGB SR stars (Figs.[Tl), because in the spectra of 
Miras there is strong water absorption. RSGs and AGB SR stars 
have w eaker water absorption than AGB Miras. lYang & Jiand 
(l201lh suggest that RSGs have bluer [3.6] - [4.5] due to a con- 
tinuum depression around 4.5 //m by CO bands. RSG-SR stars 
have redder [3.6] - [4.5] colors than RSGs with irregular pulsa- 



tions. Only weak water absorption is seen in some variable RSG 
stars. 

In Fig. [3 we plot our synthetic GLIMPSE colors of RSG 
stars versus mass-loss rates, as inferred by IVerhoelst et al.l 
(I2OO9I) . Both [3.6] - [4.5] and [3.6] - [8.0] colors correlate with 
mass-loss rates. The correlation of mass loss with the [3. 6] -[8.0] 
is explained by the increasing silicate strength at 9.7 fim. The 
correlation of mass loss with [3.6] - [4.5] color can be due to 
CO at 4.3 jum, and/or to the existence of an extra source of op ac- 
ity (dust component), as suggested bv IVerhoelst et alj (l2009h . 
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Fig. .1. Left: A representative samples of ISO-SWS spectra of AGB stars. Over-plotted transmission curves of the four 
SPITZER/IRAC filters. Right: A representative samples of ISO-SWS spectra of RSG stars. Over-plotted transmission curves of 
the four SPITZER/IRAC filters. 
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